Route to Lambda in conformally coupled phantom cosmology 
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In this letter we investigate acceleration in the flat cosmological model with a conformally coupled 
phantom field and we show that acceleration is its generic feature. We reduce the dynamics of the 
model to a 3-dimensional dynamical system and analyze it on a invariant 2-dimensional submanifold. 
Then the concordance FRW model with the cosmological constant A is a global attractor situated 
on a 2-dimensional invariant space. We also study the behaviour near this attractor, which can 
be approximated by the dynamics of the linearized part of the system. We demonstrate that 
trajectories of the conformally coupled phantom scalar field with a simple quadratic potential crosses 
the cosmological constant barrier infinitely many times in the phase space. The universal behaviour 
of the scalar field and its potential is also calculated. We conclude that the phantom scalar field 
conformally coupled to gravity gives a natural dynamical mechanism of concentration of the equation 
of state coefficient around the magical value w c s — — 1 . We demonstrate route to Lambda through 
the infinite times crossing the w e g = — I phantom divide. 

PACS numbers: 98.80.Bp, 98.80.Cq, 11. 15. Ex 



At present the scalar fields play a crucial role in modern cosmology. In an inflationary scenario they generate an 
exponential rate of evolution of the universe as well as density fluctuations due to vacuum energy. The Lagrangian 
i— i , for a phantom scalar field on the background of the Friedmann-Robertson- Walker (FRW) universe is assumed in the 
form 

m ■ 

in ' = llg^d^W + W - 2EW)], (f) 

m : 

where is the metric of the spacetime manifold, ip = ip(t), t is the cosmological time, R = R{g) is the Ricci scalar 
for the spacetime metric*?, £ is a coupling constant which assumes zero for a scalar field minimally coupled to gravity 
■^J" ' and 1/6 for a conformally coupled scalar field, U(ip) is a potential of the scalar field. 

The minimally coupled slowly evolving scalar fields with a potential function U(ip) are good candidates for a 
description of dark energy. In this model, called quintessence [l|, [|| , the energy density and pressure from the scalar 
field are = — l/2vj 2 + U(ip), p^ — —l/2ip 2 — U(ip). From recent studies of observational constraints we obtain that 
Wip = p^/p^, < —0.55 [3|. This model has been also extended to the case of a complex scalar field [J, HI]. 

Observations of distant supernovae support the cosmological constant term which corresponds to the case ijj ~ 0. 
Then we obtain that = — 1. But there emerge two problems in this context. Namely, the fine tuning and the cosmic 
coincidence problems. The first problem comes from the quantum field theory where the vacuum expectation value 
is of 123 orders of magnitude larger than the observed value of 10~ 47 GeV 4 . The lack of a fundamental mechanism 
which sets the cosmological constant almost zero is called the cosmological constant problem. The second problem 
called "cosmic conundrum" is a question why the energy densities of both dark energy and dark matter are nearly 
equal at the present epoch. 

One of the solutions to this problem offers the idea of quintessence, which is a version of the time varying cosmological 
constant conception. Quintessence solves the first problem through the decaying A term from the beginning of the 
Universe to a small value observed at the present epoch. Also the ratio of energy density of this field to the matter 
density increases slowly during the expansion of the Universe because the specific feature of this model is the variation 
of the coefficient of the equation of state with respect to time. The quintessence models @, @ describe the dark energy 
with the time varying equation of state for which wx > — 1, but recently quintessence models have been extended to 
the phantom quintessence models with wx < — 1- In this class of models the weak energy condition is violated and 
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such a theoretical possibility is realized by a scalar field with a switched sign in the kinetic term ip 2 — » — tp 2 0, H, [I]- 
From theoretical point of view it is necessary to explore different evolutional scenarios for dark energy which provide 
a simple and natural transition to Wx = — 1. The methods of dynamical systems with notion of attractor (a limit 
set with an open inset) offers the possibility of description of transition trajectories to the regime with wx = — 1- 
Moreover they demonstrate whether this mechanism is generic. 

Inflation and quintessence with non- minimal coup ling constant are studied in the context of formulation of necessary 
conditions for the acceleration of the universe jlOl ] (see also [TJ E3)- We can find two important arguments which 
favour the choice of conformal coupling over £ ^ 1/6. The first, equation for the massless scalar field is conformally 
invariant [HI [3. The second argument is that if the scalar field satisfy Klein- Gordon equation in the curved space 
then tp does not violate the equivalence principle, and £ is forced to assume the value 1/6 [l5[ . 

While recent astronomical observations give support that the equation of state parameter for dark energy is close to 
constant value — I they do not give a "corridor" around this value. Moreover, Alam et al. [l6| pointed out that evolving 
state parameter is favoured over constant wx = — 1 • The first step in the direction of description of the dynamics of 
the dark energy seems to be investigation of the system with evolving dark energy in the close neighbourhood of the 
value wx = — 1. For this aim we linearize dynamical system at this critical point and then describe the system in a 
good approximation (following the Hartman-Grobman theore m 1171) by its linearized part. 

Other dark energy models like the Chaplygin gas model [TltiaH, @, and references therein] and the model with 
tachyonic matter can also be interpreted in terms of a scalar field with some form of a potential function. 

Recent applications of the Bayesian framework [2l|, |22|, |23j, |24j, |25| of model selection to the broad class of cosmo- 
logical models with acceleration indicate that a posteriori probability for the ACDM model is 96%. Therefore the 
explanation why the current universe is such close to the ACDM model seems to be a major challenge for modern 
theoretical cosmology. 

In this letter we present the simplest mechanism of concentration around wx = — 1 basing on the influence of a 
single scalar field conformally coupled to gravity acting in the radiation epoch. Phantom cosmolo gy non-minimally 
coupled to the Ricci scalar was explored in the context of superquintesence (wx < —1) by Faraoni [H, H3] and there 
was pointed out that the superacceleration regime can be achieved by the conformally coupled scalar field in contrast 
to the minimally coupled scalar field. 

Let us consider the flat FRW model which contains a negative kinetic scalar field conformally coupled to gravity 
(£ = 1/6) (phantom) with the potential function U(ip). For the simplicity of presentation we assume U(ip) oc tp 2 . 
In this model the phantom scalar field is coupled to gravity via the term ^Rtp 2 . We consider massive scalar fields 
(for recent discussion of cosmological implications of massive and massless scalar fields see [ll])- The dynamics of a 
non-minimally coupled scalar field for some self-interacting potential U(tp) and for an arbitrary £ is equivalent to the 
action of the phantom scalar field (which behaves like a perfect fluid) with energy density p^ and pressure [5^] 

H = ~\^ + U W - 3 £#V ~ 3£ff(V 2 ); (2) 



1 



v?-c/(v)+e 2ff(v 2 y+(v 2 )' 



2H + ZH Z 



(3) 



where the conservation condition p^ = — 3-ff (p,/, -I- p^) gives rise to the equation of motion for the field 

tp + 3Hip + £Rip 2 -U'(ip) =0, (4) 

where R = 6(ij + 2i7 2 ) is the Ricci scalar. 

Let us assume that both the homogeneous scalar field tp(t) and the potential U(ip) depend on time through the 
scale factor, i.e. 



^(t)=^(a(i)), U(ip) = U(iP(a)); 



(5) 



then due to this simplified assumption the coefficient of the equation of state w^ is parameterized by the scale factor 
only 



and 



w^ = W-,p(a), p,p = w^(a)p^(a), 
4>' 2 H 2 a 2 - U(i/>) +£[2(iJj 2 yH 2 a + (^ 2 )] + £ [H + 3i? 2 ]^ 2 



\^)' 2 E 2 a 2 + U(ip) - 3£H 2 tP 2 - 3^ 2 )'H 2 c 



(6) 



(7) 
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where prime denotes the differentiation with respect to the scale factor. 

We assume the fiat model with the FRW geometry, i.e., the line element has the form 



ds 2 = -dr 



l (d8 2 + sin 2 6dip 2 )], 



(8) 



where < <p < 2ir, < 6 < tt and < r < oo are comoving coordinates, t stands for the cosmological time. It is also 
assumed that a source of gravity is the phantom scalar field ip with the conformal coupling to gravity £ = 1/6. The 
dynamics is governed by the action 



s = \ /<i>, v — 



1 



m 2 p R + {sT^ v + -RV - 2UW)) 



(9) 



where m 2 = (8irG) 1 ; for simplicity and without lost of generality we assume 47rG/3 = 1 and U(ip) is the scalar field 
potential 



U{i/>) = -m 2 ^ 2 . 



(10) 



After dropping the full derivatives with respect to time, rescaling phantom field if) — * (f> = tpa and the time variable 
to the conformal time dt = adr/ we obtain the energy conservation condition 



c 1 '2 i 1 J./2 1 2 2/2 

I = -a + -<p - -m a <p = p rfi 



where p r $ is constant corresponding to the radiation in the model. The equations of motion are 



(11) 



4>" = m 2 a 4> 



(12) 



where a prime denotes the differentiation with respect to the conformal time dt = adr] and m 2 > 0. 
From the energy conservation condition we have 



1 



1™2„2 ±2 



Pr,0 + 2 771 a 
1 



(13) 



and now from the equations of motion (|12| we receive 



(prfi + -rn 2 a 2 <f) 2 )(j) + -m 2 a<f)(l + 4> 2 



The effective equation of state parameter is 



for our model this parameter reduces to 



Weff 



P4> + zPr 
P<P + pr 



a) = 0. 



Wcff 



1 \(j}' 2 + \m 2 a 2 cj) 2 - p rfi 
3 -\$ 2 + hn 2 a 2 (j) 2 + p rfi 



(14) 



(15) 



(16) 



where a prime denotes the differentiation with respect to the conformal time and finally taking into account equa- 
tion (fT3"|) we have 



W c f£ 



+ \m 2 a 2 ^) 2 - prfi 
^m 2 a 2 cj) 2 + p r:0 



(1 + ^ 2 ) 



(17) 



For a,4>^> prfi this equation reduces to 



™ cff = --(20 2 + l), 



(18) 



and it is clear that for any value of <f> w e s is always negative. 
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To analyze equation (|14p we reintroduce the original phantom field variable tp = — and da /a — din a. Now equation 
(HU) reads 



(19) 



where a prime now denotes the differentiation with respect to a natural logarithm of the scale factor. Introducing 
new variables y = ip' and p r — p r fla~ A we can represent this equation as an autonomous dynamical system 



4>' = y 



p' r = -4p r 



hm 2 ip 



\m 2 ip 2 



(V(y + V-)-i)(i + (2/ + ^) 2 ) 



(20) 



There are the two critical points in the phase space (ip,y,p r ), namely tp — ±1, y — 0, p r = 0. The linearization 
matrix reads 



.4 : 



1 

-2(1 + ^ 2 ) -1-(1+^ 2 ) ; -| F (V, 2 -l)(l + ^) 



V=0, Pr =0 



1 

-4 -3 




(21) 



The eigenvalues for this matrix are Ai,2 = §(— 3 ± «v7) and A3 = —4. 

To find a global phase portrait it is necessary to study the system in the neighbourhood of the critical points which 
correspond, from the physical point of view, stationary states (or asymptotic solutions). Then the Hartman-Grobman 
theorem guaranties us that the linearized system at this point is a well approximation of the nonlinear system. First, 
we must note that p r = is in the invariant submanifold of the 3-dimensional nonlinear system. It is also useful to 
calculate the eigenvectors for any eigenvalue. We obtain following eigenvectors 



Vl,2 



3 " 

8 




r ^ ' 




"(f 


1 




8 




, ^3 = 















1 



(22) 



They are helpful in construction of the exact solution of the linearized system 

3" 4 * 
V7 
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0" 
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-3 
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-4 




1 
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e 2 cos-^-t — e 2 sin 
e~5*sin-^i e^'cos-^-i 



2 L 



1 

1 



_8_ _3_ 



XQ 

yo 



(23) 

where x = tp — ipo, y — tp' — ip' , z = p r — po and xq, yo, zq are initial conditions and we have substituted In a = t. 

If we consider linearized system on the invariant stable submanifold z = 0, it is easy to find the exact solution. If 
we return to the original variables tp,tp', then ip(\na) is the solution of the linear equation 



(*p - iM" + 3{ip - tp y + 4(v> - iM = 0, 



(24) 



( (. ■ — < 1 ) = ( ' I < wp ( — — In a \ cos ( In a J + C'2 cxp f— — In a \ sin ( In a 



(25) 



i>o) = C\a 2 cos i In a j + C^a 2 sin In a 



(26) 



Because of the lack of alternatives to the mysterious cosmological constant [2l|, |22j we allow that en ergy might vary 
in time following assumed a priori parameterization of w(z). In the popular parameterization (30l. l3ll. 1321 133| appears 
a free function in most scenarios which is a source of difficulties in constraining parameters by observations. However 
most parameterizations of the dark energy equation of state cannot reflect real dynamics of cosmological models with 
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dark energy. The assumed form of w(z) can be incompatible with the w{z) obtained from the underlying dynamics 
of the cosmological model. For example some of parameters can be determined from the dynamics which can be 
crucial in testing and selection of cosmological models [21]. Our point of view is to obtain the form of w(z) specific 
for given class of cosmological models from dynamics of this models and apply it in further analysis both theoretical 
and empirical. In practice we put the cosmological model in the form of the dynamical system and linearize it around 
the neighbourhood of the present epoch to find the exact formula of w{z). For the phantom scalar field model this 
incompability manifests by the presence of a focus type critical point (therefore damping oscillations) in the phase 
space rather than a stable node (Fig. [1] and its 3D version Fig. [2]). 

The properties of the minimally coupled phantom field in the FRW cosmology using the phase portrait have been 
investigated by Singh et al. [HI (see also [35[ for more recent studies). Authors showed the existence of the deSitter 
attractor and damped oscillations (the ratio of the kinetic to the potential energy \T/U\ to oscillate to zero). 

We can also express w e g in these new variables 

=~hi>+ v/) 2 - Pr ~ Yf l2 m } (27) 

o { p r + im^V J 



dw off 2 f m 2 V> 2 /i 1 1\ / 1 1 t i tt\ i 2,2 2 + 



din a 3 |_ p r + ^m 2 ip 2 



(yj + W + r) + fr-mV , ^tt/2^2 I 1 + U> + ^') 2 ] } • ( 28 ) 

2'"- y (Pr + 2 m r J ) 

Recently Caldwell and Linder [3 6] have discussed dynamics of quintessence models of dark energy in terms of w — w' 
phase variables, where w' was the differentiation with respect to the logarithm of the scale factor. These methods 
were extended to the phantom and quintom models of dark energy [37l l3l| . Guo et al. dU examined the two-field 
quintom models as the illustration of the simplest model of transition across the wx = — 1 barrier. The interesting 
mechanism of acceleration with a periodic crossing of the w = — 1 barrier have been recently discussed in the context 
of the cubic superstring field theory [39J. In the model under consideration we obtain this effect but trajectories cross 
the barrier infinitely many times. The main advantage of the discovered road to A is that it takes place in the simple 
flat FRW model with the quadratic potential of the scalar field. 

It is easy to check that at the critical points w c s = — 1 and grjff = 0. Since these points are sinks there is infinite 
many crossings of w c g = — 1 during the evolution. 

The methods of the Lyapunov function are useful in discussion of stability of the critical point of the non-linear 
system. The stability of any hyperbolic critical point of dynamical system is determined by the signs of the real parts 
of the eigenvalues Xi of the Jacobi matrix. A hyperbolic critical point is asymptotically stable iff real A.; < Vz, if xq 
is a sink. The hyperbolic critical point is unstable iff it is either a source or a saddle. The method of the Lyapunov 
function is especially useful in deciding the stability of a non-hyperbolic critical points [TJ, p.129]. The construction 
of the Lyapunov function was used by [4p| for demonstration that periodic behaviour of a single scalar field is not 
possible for minimally coupled phantom scalar field (see also [4l|). 

The quantity w' cS in terms of w e ff and (ln-0)' reads 

w' eS = -(1 - 3w ff)(l + Wcff + ^(ln^)'). (29) 

It is interesting that equation (|29p can be solved in terms of w(a) - the mean of the equation of the state parameter 
in the logarithmic scale defined by Rahvar and Movahed [13] as 

/>( ffl ')d(W) 
Ji d(lna') 



namely: 



>(a) = I_i a 3(i+*(°V. (31) 



They argued that this phenomenological parameterization removes the fine tuning of dark energy and px/Pm oc 
a -3w(a) approaches a unity at the early universe. Note that in w(a) = — 1 that 



4) + i = J(i-A 



(32) 




FIG. 1: The phase portrait (w e ff, w' e g) of the investigated model on the submanifold p r = 0. This figure illustrates the evolution 
of the dark energy equation of the state parameter as a function of redshift for different initial conditions. In all cases trajectories 
cross the boundary line w e « = —I infinite many times but this state also represents the global attractor. 



where 



. = < u - ( I H- : cos(^ ln(l + z))-C 2 sin(^ ln(l + z)) 



(33) 



In Fig. [3] we present the relation w(z) for different values of parameters ipo = ±1, C\ and Ci. 

In this letter we regarded the phantom scalar field conformally coupled to gravity in the context of the problem of 
acceleration of the Universe. We applied the methods of dynamical systems and the Hartman-Grobman theorem to find 
universal behaviour at the late times - damping oscillations around u> e g = — 1. We argued that most parameterizations 
of the dark energy, such as linear evolution of w(z) in redshift or the scale factor, cannot reflect realistic physical 
models because of the presence of non- hyperbolic critical point of a focus type on the phase plane (w,w'). We 
suggested a parameterization of a type 



w x (z) = -1 + (l + z) 3 \ CiCOs(ln(l + z)) + C 2 sin(ln(l + z)) 



(34) 



which parameterizes damping oscillations around wx = —1 "phantom divide", and finally, with the help of this 
formula one can simply calculate energy density for dark energy px 



Px = Px,oexp ( - B) cxp ((1 + z) 3 [^sin(ln(l + z)) + flcos(ln(l + z))]). 



(35) 
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